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We investigate the influence of the R-parity violating couplings A, A' and A" on the branching 
ratio of b — > sj in leading logarithmic approximation. The operator basis is enlarged and the 
corresponding 7-matrix calculated. The matching conditions receive new contributions from the 
R-parity violating sector. The comparison with the experiment is rather difficult due to the model 
dependence of the result. 
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I. INTRODUCTION 



The decay b — ► S7, forbidden at tree-level in the Standard Model (SM), is an excellent candidate for explor- 
ing the influence of new physics beyond the SM. However, the experimentally measured branching ratio Br(£> — > 
X sl ) = (3.15 ± 0.93) • 1(T 4 § is in perfect agreement with the SM prediction computed at the next-to-leading 
■ order:Br(6 — > sj)sm — (3.28 ± 0.30) ■ 10~ 4 0. This leads to the conclusion that the influence of new physics on this 
decay is either very limited or the new contributions cancel among each other to a large extent. The most serious and 
attractive extention of the Standard Model is Supersymmetry (SUSY). It has a variety of very appealing features. 
' For instance it provides a natural solution to the hierarchy problem. SUSY doubles the particle spectrum of the SM 
providing every fermion with a bosonic partner and vice versa. If SUSY were an exact symmetry of nature the new 
particles would be of the same mass as their partners. This is definitely excluded from experiment. Therefore SUSY 
must be broken. To retain the solution of the hierarchy problem one allows only a breaking which does not introduce 
quadratic divergences in loop diagrams. Although this reduces the set of possible breaking terms substantially to the 
so-called soft-breaking terms the number of free parameters of softly broken SUSY still exceeds one hundred. This 
theory, i.e., SUSY with soft breaking terms and no others, emerges naturally as low energy limit of local supersymmc- 
1 try, Supergravity (SUGRA) by breaking it at some high scale ~ 10 11 GeV and taking the so-called flat limit |§. The 
connection to supergravity eliminates much of the freedom in choosing the parameters for the soft breaking terms, 
1 enhancing the predictive power of the model. 

Viewing naturalness as a first principle in model building one runs immediately into a problem: The Yukawa interac- 
• tions are fixed by the so-called Superpotential W, which is the most general third degree polynomial that can be built 
of gauge invariant combinations of the (left-handed) superfields of the theory. In the case of the minimal extension of 
the standard model with gauge group SU(3) x SU(2) x U(l) it shows the form 



w = XfcQiHiU? + x!jiQ,ir, + a;,//,/., a; + [iH x h 2 

+mUH2 + ^a,,/. /.,/., /•.;: + x^.L.Qjr, + ±M? jk u?Dp c k , (i) 

where Q, U c , D c , L, E c , Hi and H2 label the left-handed superfields that describe left- and right-handed (s)quarks 
and (s)leptons and the Higgs-bosons (-fermions) respectively. The terms on the second line of (|l|) lead to unwanted 
baryon- and lepton- number violating vertices. The product A^- fc • A^- fc , for instance, is restricted to be smaller than 
10~ 10 (I- It is not obvious why the coefficients of these terms should not be of order unity if the corresponding 
interaction is not protected by any symmetry. Hence one way to avoid these unwanted couplings is the invention of 
a new discrete symmetry, called R-Parity j^]. The multiplicative quantum number R is then defined as 

R=(-l) 3B+L+2s , (2) 

where s is the spin of the particle. The particles of the standard model are then R-parity even fields while their 
supersymmetric partners are R-parity odd fields. The superfields adopt the value of R from their scalar components. 
The terms on the second line of (0) are then forbidden by this symmetry. One ends up with the Minimal Super- 
symmetric Standard Model (MSSM) §-§. One could even go one step further and promote this new symmetry to a 
U(l) gauge symmetry (R-symmetry). The R-charges must then be chosen such that the "nice" terms remain in the 
Lagrangian and the unwanted terms will not be allowed anymore. Further requirement is the vanishing of possible 



1 



additional anomalies [p|JlO|l . This scenario is preferred from a string theoretical view because in string theory one has 
many additional t/(l)'s "floating around", what makes the introduction of this symmetry more natural. 
It is interesting to explore what the constraints on the R-parity violating couplings, especially Xijk, Kjk> anc ^ ^ijk< are 
from the experimental point of view. Bounds on these couplings have been found by many authors [|IT|-|l4|] . Sometimes 
the bounds on products of couplings are more restrictive than the products of the individual bounds. Depending on 
the reactions that the couplings are involved in the constraints from experiment are very strong or rather poor. 
In this paper we want to explore the (theoretical) influence of A, A', and A" on the decay b — > sj. This has been 
done before by other authors Jl5|. However, here we will include the full operator basis at leading-log of the effective 
low-energy theory. 

The comparison with the experiment must result in bounds that arc model dependent. The relevant couplings are 
proportional to the inverse mass-squared of particles (Higgs, SUSY-partners) that have not yet been detected. How- 
ever, it is clear from next-to-leading-log calculations that the prediction of the MSSM with a realistic mass spectrum 
lies within the current experimental bounds. 

Because of the strong model dependence we do not perform our calculations at highest precision. Nevertheless we try 
to include all the possibly relevant contributions at leading-log. We do not claim our results to be very accurate. It 
is the aim of this paper to explore if b — ► sj has the potential power to reduce the bounds of (products of) some of 
the R-parity breaking couplings substantially, i.e., by some order of magnitude. 

This article is divided as follows: Section 2 introduces the model we are working with. We try to describe as precise 
as possible what our assumptions are. The following section deals with the effective Hamiltonian approach. The 
enhanced operator basis is presented and the 7-matrix as well as the matching conditions at Mw calculated. The 
comparison with the experiment is performed in section 4. Section 5 contains our conclusions. In the appendix we 
present some technical details of our computations, namely the mixing matrices, the relevant part of the interaction 
Lagrangian and the RGE's that are needed. 

II. FRAMEWORK 

In supersymmetry the matter fields are described by left-handed chiral superfields E*. They contain a scalar 
boson z % and a two-component fermion ip l . Real vector superfields V a are needed to form the gauge bosons 
and the gauginos A a . The minimal supersymmetric standard model is the model with the smallest particle content 
that is able to mimic the features of the standard model, i.e., including all the observed particles, gauge group 
5C/(3) co i our x 5t/(2) wca k x U(1)y> spontaneous symmetry breaking and the Higgs mechanism. Its superfields together 
with their components are collected in table [|. A few comments are in order: 

• Because the theory only deals with left-handed chiral superfields the SU (2)-singlet matter fields must be defined 
via their charged conjugated (anti-) fields. 

• No conjugated superfields are allowed in the superpotential W. Therefore we need to introduce a second Higgs 
field to give the up and the down quarks a mass when the neutral components gain a vacuum expectation value 
(vev). 

• The L's and R's in the names of squarks and sleptons only identify the fermionic partners. These fields are just 
normal complex scalar bosons. 

It must be mentioned that the fields in table | are not the physical fields. 

• In the Higgs sector three degrees of freedom are eaten by the gauge bosons in analogy to the SM. We end up 
with one charged and three neutral Higgs bosons. 

• Higgsinos and gauginos of SU(2) x U(l) mix to form charginos and neutralinos. 

• Photon, W- and Z-boson form when the electroweak symmetry breaks down. 

• The three generations of quarks and leptons mix via the Cabibbo-Kobayashi-Maskawa-matrix K to give the 
mass eigenstates in complete analogy to the standard model. 

• The family mixing also takes place in the squark and slepton sector. However, there is an additional mixing 
between the partners of left- and right-handed fermions due to the soft-breaking terms. 
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Appendix ^ gives a more detailed description of the different mixings including a complete listing of all relevant 
mixing matrices. 

The component expression of the Lagrangian we base our model on can be written as 

£ = £kin + Ant + £\V — Koft , (3) 

where 

• £kin and C m t stand for the kinetic energies, the interactions between chiral and gauge fields and part of the 
scalar potential. 

• Cw contains the rest of the scalar potential and the Yukawa interactions: 
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, 9 W i t. 



(4) 



where W should be viewed as a function of the scalar fields. 

• W is the superpotential which contains all possible gauge invariant combinations of the left-handed superficlds 
(but not their conjugated right-handed partners). In our case, this results in equation (|l|). It should be 
noted that every term is a gauge invariant combination of the corresponding superfields, for instance, QiU^U^ 

is abbreviated for e a p Qf A U^ A , where a and /3 are S'i7(2)-indices, A is an SU (3) -index and e a p is the 
completely antisymmetric tensor with 6x2 = 1 | The second line of (Q) represents the R-parity breaking 
sector. As a consequence of the antisymmetry in the fields A is antisymmetric in its first two indices and 
A" is antisymmetric in its last two indices. Therefore the trilinear R-parity breaking couplings of W contain 
9 + 27 + 9 = 45 new parameters. 

• V^oft includes the soft-breaking trilinear terms of the scalar potential and mass terms for the scalar fields and 
the gauginos. It has the form 



+jxhji2 + ijL2iiih 2 ] + h.c. + m 2 ab z a z bji 
1 



M [A a A a +h.c.]. (5) 



The MSSM in its full generality with an additional R-parity breaking sector involves over 150 free parameters. 
These are by far too many for the model to be predictive. In the following we will reduce the parameter space 
substantially by making some assumptions that are, hopefully, well motivated. 

As a first step it is important to mention that we see our softly broken global SUSY at a low energy scale ~ Mz 
emerging from a spontaneously broken local supersymmetry at a high scale Mx ~ 10 16 GeV taking the flat limit 
-Mpianck — > oo, ttiq := m grav iti n o ^constant. This fixes most of the parameters of V so ft at Mx- 

• All the coefficients of the trilinear terms in V so ft are related to the corresponding terms of W by a multiplication 
with a universal factor Am^: 



x The antisymmetry of the SU(2) product is the reason why a term ~ H\H\E1 is not introduced 
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• An analogous statement holds for the bilinear terms: 

p, = Bm fi fan = Bm (j,2i, 

where usually 

B = A-l. 



(7) 



(8) 



• The mass term of the scalars are diagonal and universally equal to the gravitino mass mo: 

m lb = mlSab 



(9) 



We assume unification of the gauge group at M x ■ As a consequence all the gaugino masses are equal at that scale: 

Mi(Mx) = M Vi (10) 

Not all entries of the Yukawa-matrices A^-, \fp and \\j are observable in the SM. One usually chooses two of them 
(in most cases X d and A e ) to be diagonal. Although this is in principle not possible in our model we will adopt this 
choice here for convenience. All the entries at Myy are then fixed by the quark/lepton masses, the vevs v% and v 2 of 
the neutral Higgs bosons Hi and H% respectively and the Cabibbo-Kobayashi-Maskawa-matrix K. 

[iH\H2 is the so-called /^-term. The mass parameter fi must be of order of the weak scale whereas the natural 
scale would be the Planck mass Mp ~ 10 19 GeV. The question why this parameter is so small is referred to as the 
fi-problem. 

/Z2iiji?2 and fliitihi mix Higgs and leptonic sector. We choose to set \i-u{Mx) = &2i{Mx) = 0. At Mz, ^n^iHi 
can be rotated away with the help of a field redefinition of the Higgs field whereas Ji-a is, at least in the case of a 
physically realistic spectrum, small enough to be neglected. 
One ends up with the following free parameters: 



A, m , M, fj, 



(11) 



Usually, one replaces one of these parameters by tan (3 = V2 jv\ . A second parameter will be fixed by the requirement 
of a correct electroweak symmetry breaking. This means, the minimum of the scalar Higgs potential must occur at 
values (vi, 1*2) which reproduce the correct mass of the Z-boson: 



Ml 



9l 



9l 



(vi + 4) 



It has been realized by many authors |l6|-|l§|] that the tree-level potential 



V a = (^ + m 2 HlHl )(h1) 2 + {y 2 + m 2 H2H2 )C h 2) 2 + 2M? h° 2 



(12) 



(13) 



9l + 9l 



is not enough to gain sensible values for v\ and v 2 - Thus we have to include the first correction AV to the effective 
potential. At a mass scale Q, it has the form |p^ , po| 



AV = 



1 

64tt 2 
1 

64tt 2 



Str 



M 4 In 



Q 2 



^(-l) 2 M2 Sp +lKM p 4 [ln 

p V 



Q 2 



(14) 
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Here, Str denotes the supertrace and M 2 is the tree-level mass matrix squared, p runs over all particles of the theory 
with spin s p whereas M p is the corresponding eigenvalue (mass) of the particle. n p counts for the degrees of freedom 
according to colour and helicity. The eigenvalues M p depend on the neutral components of hi and hi and therefore 
change the shape of the potential and its minimum. For most of the particles, they can only be computed numerically. 
One comment is in order: A phase rotation of the Higgs fields can turn a negative vev into a positive one. This freedom 
is reflected in the fact that the sign of /< can be chosen freely giving then a different phenomenology. 
To compute the mixing matrices and the one-loop effective potential one has to know the mass matrices at Mz- 
Unfortunately, for most of the parameters we know the boundary conditions at the high scale Mj. Hence one has 
to set up the complete set of RGE's to run the parameters from Mx to Mz- The details of how to get a consistent 



parameter space are described in section IV, the complete set of RGE's can be found in appendix 



III. THE EFFECTIVE HAMILTONIAN 
A. The case of the MSSM 

The decay b — > 57 occurs at energies of a few GeV~ mj,. This is much below the weak scale. It makes sense to 
work with an effective Hamiltonian 7i e ff where all the heavy fields (compared to mj) are integrated out In the 
SM these are the Z- and the VT-boson and the top quark whereas for our purposes the Z-boson does not play any 
role. A result of integrating out these fields is the appearance of new local operators of dimension higher than four. 
This can be illustrated by the shrinking of the Feynman diagrams in Fig. |l|. 

Once one lets the strong interaction come into the game, QCD corrections of these new operators give rise to additional 
operators. For a consistent treatment of all the corrections at a certain level of QCD we have to include a set of 
operators Oi which closes under these corrections. Our low-energy theory is then described by an effective Hamiltonian 

H eS = -J2 c i°i- ( 15 ) 

i 

The QCD renormalization of the operators must be performed at a scale where no large logarithms appear, i.e., 
at ~ My/. However, a calculation of Br(6 — * sj) at energy scales ~ nib requires the knowledge of the Wilson 
coefficients Cj at that scale. The Cj's depend on a renormalization scale [i. They obey the renormalization group 
equations 

' ICi =1*03, (16) 



din /i 



where 7^ is the gamma-matrix that emerges from the QCD-mixing of the operators Oi- As initial conditions we find 
Ci(My/) by matching the effective theory with the full theory at that scale. The RGE's are then solved to give Ci{rrib) 

In the standard model the relevant set of operators for b — > sj is given by 

0\ = (sTa^bLa) (cl pl^CL/l) 

02 = (slaJ^bLfi) (cZplnCLa) 

03 = (~sZ a j' 1 b La ) 22(qLii37»qLi0) 

i 
i 
i 
i 

O7 = ——^ m b ~sl a a^ v b Ra F^ v 

WTT Z 
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° 8 = if? mb ^°^b R0 t aa ? G a ^. (17) 

Here, Ql/rI^Il/r = 97^(1 T 7 5 )<7, SLC^b^ — say,„(l + 7 5 )6, the sum runs over the five quarks of the effective 
theory at mj,, and a, /? are colour indices. The operators 0\ and O3 — Og are introduced by QCD corrections through 
diagrams like those depicted in Fig. ||. 

It is interesting to note that in the MSSM without R-parity violation the relevant operator basis does not change 
although we definitely have new decay channels. The 7-matrix at leading-log can be found in appendix ^ by picking 
up the relevant entries. Because the contribution of the diagram in Fig. ^| is not divergent, the mixing of 0\ — Oq 
with O-j /§ involves two-loop diagrams like those of figures |J and || |^,^3[ . 

Although we only need the divergent part of these diagrams one has to be very careful in computing the countert- 
erms because one must include certain additional operators pi| , ^5f , so-called evanescent operators that vanish in four 
dimensions but must be kept in D dimensions in all intermediate steps of the calculation. Furthermore, the two-loop 
results are regularization scheme dependent pEfl . This regularization scheme dependence of the 7-matrix cancels with 
possible finite one-loop [but 0(a®)] contributions from O5 and Oq, inserted in the diagram of Fig. |^, to the matrix 
element of b — > sj. As a result of these complications, the calculation of the 7-matrix at leading- log has been finished 
only few years ago p6|-|29|] . Meanwhile, also the next-to leading result is known pC^- |43| , p5| ] . We will concentrate on 
the leading-log calculation. 

The matching of the effective theory with the full theory at My/ must be performed only at order a° because the 
leading QCD corrections are already included in the operator mixing. In the standard model this involves diagrams 
with a VK-exchange for O2 and diagrams with a W-t-loop (in the unitary gauge) for 7 / 8 corresponding to Fig. || (a). 
The result is then Q 

C 2SM (M W ) = -^K* ts K tb (18) 

C 7SM (M W ) = -^KlKtbSxtwiQuF^xtw) + F 2 {x tw )} (19) 
G F 

Cssm(M w ) = --j=K ts K tb 3x t wFi(x tw ), (20) 

where Gf = c&/{Ay/2My V ) w 1.166 • 10~ 5 GeV -2 is the Fermi constant, x a b = m^/rn 2 and the functions Fi are given 
in appendix y. 

The matching conditions for C?/g become far more complicated in the case of the MSSM, even without R-parity 
violation. In addition to the W-t-loop there are four more combinations of particles in the loop: 

• charged Higgs H ± — top t ,Fig. ^| (b) 

• up-squark u — chargino x ch ,Fig- H (c) 

• down-squark d — neutralino x° jFig- 1§ (d) 

• down-squark d — gluino g ,Fig. ^ (e) 

In principle, these contributions give rise to two new operators 

~ 6 

7 = ^T"2 m b ~SRa°~nvb La F^ V 

Note that O7/8 differ from only by their handedness. The Wilson coefficients of these new operators are 

usually so small that they can be safely neglected ]4q] . However, we will include them in our calculations because we 
need them later on anyway. We only neglect contributions from the Higgs sector which are proportional to some light 
quark masses. 

The matchi ng co nditions become rather involved now. They include many mixing matrices, whose definitions we give 



in appendix A 1 . We have 
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CVmSSM = C7SM 



- ^K; s K tb {cot 2 [3 x tH [QuFi (x tH ) + F 2 (x tH )} 
+ x t H[QuF 3 (x tH ) + F 4 (x tH )}} 



I Dei ; > 

AmT 23 1 3ji 

Uj 

1 1 m YS h 



ryd Ad* 



Fiix^.) + Q u F 2 {x x c^.) 



Qd 1 

4 m 2 - 

a a 



D% t D% lF ^x x% ) + -^D^Cl^ix^] 



2 ^ 2 1 

3 m~ 



(22) 



CVMSSM = CqsM 



--^K t * s K tb [cot 2 (3 xt H F 1 {x tH ) + x tH F 3 (x tH )} 
v2 



1 1 



+ 



4 «4 3 
1 1 

4m 2 ; 



B. 



2je B 3jt F 2(x x oh il .) + -^-B^A^F^x^h^.) 



Ad* 



nib 



2 2 1 

+ 3 5 W 



r%Ti j3 F 2 (x g3j ) - ^ 2j T% 3 F 4 (x gii ) 



3 2 1 

Gf 



^ 2j ri j3Fl ( Xg3j ) - Jr^t .r^Fa^- ) 



C 7 mssm = -—j=K* ts K th \ tan 2 ft xtH[QuFi(x t H) + F 2 (x tH )} 
v 2 Ti t 



-i-i— — 4 d 4 d * 



+ 



4< 

1 1 m oh 

4 m? m b A W B W 



Fi{x x f il] ) + QuF 2 {x xtil] ) 

F 3{x x c^.) + QuFi(x x c^.) 



Qd 1 

4 m 2 - 

a 7 - 



2 n 2 1 



(23) 



(24) 



Csmssm = ^Kl s K t b — s - 2 ^- tan 2 (3x t HFi(x t H) 



V2 
1 1 



A 2 jeA 3 * e F 2 (x x < : h il .) + —^-A 2je B 3 j e F 4 (x x ch il .) 



4 "4, - 
1 1 
4 m 2 - 



C d 2je Cl; e F 2 (x x0e3j ) + -^C^Df^ix^) 



2 o 1 
+ 3 5 W 



rftAi^) - ^r^ri j3 F 4 (x^.) 



7 



3 2 1 



m 



(25) 



B. Including the R-parity breaking terms 

1. The new operator basis 



As mentioned before, the difference between the standard model and the MSSM does not lie in a change of the 
operator basis but rather in different matching conditions at My/- This situation changes drastically if one includes 
the R-parity breaking sector. Now our basis has to be enlarged. To find out which are the relevant new operators we 
first write down the R-parity breaking Yukawa couplings: 



C 



Yukawa^!, 



= C\ + C\> + C\" . 



(26) 



where 



C x = -A 



ijk 

2 A u'fe 



ne ~t i r>"t~t i "n e t c 

^LitK} m e £ d Rk u L m + Tf H/ utd Rk e u + ^Rik K ] m d ^i u ^ 



(27) 



v \r.^Z.A.. T- d t IT".,,.. rd Jt 



Fji>ed Rk d Lj - T£ jt did Rk v Li - Y mk d\v c Ri d Lj 
2 



+ -^K*jk F e Lei K m jel<UL m d Rk + T^u\eLid Rk + Y R \ e K m jd e u Lm e Ri 
" d dh 

Llj a l l 

1 



Y v ti v\dLjd Rk — Y'[ e jd\i'Lid R k — Y R \ e d(dL 3 v c Ri 



1 



vjk 



A ijk 



T RH uld Rj d c Lk + Y d mj d\u Rl d c Lk + Y d mk d\u Rl d 



Rii 



d L k d Rj + T m d e d Lk U Ri + T BM d l d Li u R 



(28) 



(29) 



Here, all the fields belong to the mass basis. The colour indices have been omitted. 

The next task is to build four-quark operators out of two Yukawa couplings that contribute at 0(a s ) to b — > S7. 
The boson serves as a bridge between the fermions in analogy to the W boson in the standard model. The following 
points must be taken care of: 

• The top quark is not at our disposal in the five flavour effective theory. 

• We not only need an ingoing b and an outgoing s. The two remaining quarks must be of the same type because 
one has to be able to close the loop with these fermions. 

It is clear that C\ cannot participate because it contains no squarks and semi-leptonic operators can be neglected. 
Also, C\i does not mix with Cy . As an example, we take the first term of C\> together with his hermitian conjugate. 
The situation is depicted in Fig. [?]: 



k 2 <.m% t i l 



RULr. 



k 2 - ml 2 



7;Kb3 T Lta K mbe\u Lm b R 



(30) 
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In the first step we used the fact that the two squarks have to be of the same type and the unitarity of the CKM-matrix. 
In the second step we performed a Fierz rearrangement. This is done to get the same structure (i.e., two vectors) 
for the four-fermion operator as in the standard model case. The advantages of this rearrangement will become clear 
when calculating the 7-matrix. In the last line we put the colour indices a, (5 for clarity. 

As one can see clearly, the effect of the (unitary) squark mixing matrix T e L becomes enhanced if the masses of the 
sclcctrons are very different for the three generations. If there were a mass degeneracy they would simply give a factor 
Sia- This is a general feature in our calculations. 

The operator that appears in Eq. (|3^) is of a new type. It consists of a right-handed b- and a right-handed s-quark. 
(Actually, these are two operators, one with a pair of it-quarks and the other with two c-quarks.) 

A careful investigation results in the following set of new operators: 
From C' x one gets 



-Pi 






P-2 


= {^Ral^bRp) 


(cEplfiCLa) 


P 3 


= (sKc7 M &R/3) 


{d L pl^d La ) 


Pi 


= {sRa^blip) 




P 5 


= (SRa^bRp) 


{b L plvb La ) 


P 6 


= {sLa^bLp) 


(da ^^diia) 


Pi 


= {slal^bLfj) 


(Sfl/37M s fl") 


P, 


= {SLa^bLp) 


{bRfil^bRa) 


P, 


= (sRal^bRa) 


i 


Pw 


= (SR Q 7 M &fl/3) 


i 


Pll 


= {SRal' J 'b Ra ) 


i 


Pl2 


= (sRaJ^bap) 





i 



It is worth noting that the operators P\ - Ps emerge directly from the Lagrangian and Pg - P12 are induced through 
QCD corrections. One would expect partners of Pi - P 8 with colour structure (aa)((3(3) to be introduced by QCD. 
This does not happen at leading-log (accidently) . 

C\n leads to the following additional operators: 

Rl = (sRal^b Ra ) {URp'YpURp) 

R-2 = (sRal^bRp) (uRplftURa) 

R3 = (sRaJ^bRa) (c R ~/3'Y li CR0) 

R4 = {sRaJ^bRp) (c R ~/37f l C Ra ) 







#5 = {sRal^bRa) (^7^/3) 

-Re = (sr q 7 m &R/3) {dnp-i^dRa). (32) 
Here, all the operators i?i - R§ appear already in the tree-level effective Lagrangian. 

In principle, the effective Hamiltonian also contains semi-leptonic Operators. However, these can be neglected because, 
as a consequence of their Dirac structure, they do not contribute to the decay b — ► 57. 

2. The ^/-matrix 

The whole basis now consists of 28 operators. Their QCD-mixing is described by a 28 x 28-7-matrix. It is depicted 
in appendix There are three different blocks in this matrix that have to be treated in separate ways. 

• mixing of four-fermion operators among themselves 

This block involves the one-loop diagrams of Fig. pi There are no further complications. The mixing of 0\ - 0% 
is known already since a long time 

• mixing of O7 , Os and O7, 0$ among themselves 

These entries need not to be computed. The mixing of O7 and Os is known and the new operators mix in 
exactly the same way such that the corresponding numbers can be copied. 

• mixing of the four-fermion operators with Of, Os, O7 , Os 

This task is more difficult. In principle, one has to compute the divergent part of all the diagrams of Fig. || and 
|5] together with their counterterms and the contributions of the evanescent operators p?]]. Moreover, there are 
four different types of chiralities to be inserted: (LL) (LL) , (LL) (RR) , (RR) (LL) and (RR) (RR) . For the first 
two this calculation had to be performed for the case of the standard model. The detailed results are listed in 
. The second two types of insertions are new. Fortunately, one can deduce the divergent parts of these types 
of insertions by making the following two observations: 

— The diagrams of Fig. |] contain only one fermion line. Here, it is crucial if the two quark pairs of the 
inserted operator have the same or opposite chirality. Thus, an insertion of a (RR) (RR) leads to the same 
divergence as an insertion of a (LL)(LL) as well as the divergence of an insertion of a (RR)(LL) is the 
same as for the case of (LL)(RR). 

— One must be careful with the diagrams that contain a closed fermion loop (Fig. |J). At a first sight, one 
may deduce that the divergence does depend on the chirality of the quarks running in the loop but not 
on the chirality of the s- and the 6-quark. However, this is wrong The differ ence between a right- and a 
left-handed quark in the loop results in a term Tr(±7 5 T) where T stands for a collection of at least four 
Dirac 7-matrices and the sign corresponds to a right- or left-handed quark in the loop, respectively. The 
trace leads to an e-tensor that is contracted with a 7-matrix between the external quarks. This produces 
an additional 7 s . Now note that 7 5 Pr = +Pr, whereas 7 5 Pl = —Pl- This means, we need to change the 
chirality in both fermion pairs to end up with the same mixing. Hence, if an insertion of an operator of 
type (LL)(LL) [(LL)(RR)] gives a certain contribution to 7 / s the (RR)(RR) [(RR)(LL)] operator will 
give exactly the same contribution to O^/g. 

To summarize, the mixing of the four-fermion operators with O7/8 can be deduced completely from the results of ]26| ] 
by interchanging left-handed and right-handed projectors. 

All the previous calculations involve 7 s . It is therefore clear that the results depend on the regularization scheme. 
This dependence will be cancelled by finite one- loop (but 0(a®)) contributions of some four-fermion operators Qi to 
the Amplitude A of b — > S7 through the diagram of Fig. |l Schematically, the result is then 



2 We thank the authors of P8| for clarifying this point to us. 



46|. 
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A = C 7 (s-f\0 7 \b) u . cc + C 7 (sj\0 7 \b} tree + ^ C <2* ( S 7|Qi|&)l-loop- 



An alternative |19| is to define effective coefficients in a way that A becomes 

A = Cf{sj\0 7 \b) tiee + Cf( S1 \6 7 \b) trcc . 
This can be achieved by defining four vectors {j/i}, {z{\, {j/i}, and {z{\ through 



(33) 



(34) 



(sj\Qi\b) 



1 — loop 



= :y t (s 7 \0 7 \b) t 



(sgluon|Qi|b)i_ loop =: z j; (s gluon|0 8 |b) trco 
(s7|(3i|6)i_i 00 p =: yi(sj\d 7 \b) tiee 



(sgluon|Qi|b)i_i oop 
The effective Wilson-coefficients must be defined as 



Zi(s gluOn|6 8 |b)tre 



(35) 



Cf{n) 


■■= CM 4 


■^ViCiifi) 

i 




:= C s (n) 4 


-^ZiCiin) 

i 


Cffjj.) 


:= drill) 4 


■^mCibj.) 

i 


cf(n) 


:= Cs(ji) 4 


■^ZiCiin) 

i 









The vector 



& S (V) := {c^Cfi^Cfi^Cfi^CfM} 



(36) 



(37) 



is then regularization scheme independent. Remember that the index i runs over all four fermion operators. 

The effective coefficients obey RGE's which can be derived from the RGE's for Cfc(/x), where k labels the whole set 

of operators. They are 



Cf(M) = ^CfO*), 



In jjj 



(38) 



where 



-21 



Ijk 
Jjk 
Ijk 
Ijk 



j=l z ilji 
2_,i=l Vilji 
2-ii=l z ilji 



Vjlo 7 o 7 ~ z j7o 8 o 7 


k = Or, j 


= 1,- 


.,24 


z jlo s o 8 


k = Os, j 


= !,• 


.,24 




k = Or, J 


= 1,- 


.,24 




k = O s , j 


= !,• 


.,24 




otherwise 







(39) 



For a finite contribution of an operator inserted in the diagram of Fig. || we need 

• two pairs of fermions with different chirality, i.e., operators of the form (LL)(RR) or (RR){LL) 

• a 6-quark running in the loop. 
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This reduces the possibilities to O5, 0$, P5, P%, P%i and P12. The results for {y;}, {zi}, {jji} and {Si} are then 



-| i = 5 
-2 i = 6) P 8 Zi 
otherwise 



-f i = P n 

-2 i = ft, Pi 2 5, = 

otherwise 



2 i = 5 
otherwise 



2 i = P l: 
otherwise 



(40) 



There is one more subtlety concerning the matching at M\y: In the standard model we have (s"f\H cS \b) 
C 7 (Mw)(s^\0 7 \b) . Now, more operators contribute: 



( S1 \H ctt \b) = C 7 {M w ){ sl \0 7 \b) + C 7 {M w ){ sl \0 7 \b) 

~2C Pa (M w )( sl \P s \b) - 2Cp 5 (M w )( sl \P 5 \b) 
= Cf{M w )( sl \0 7 \b) + Cf{M w )( sl \6 7 \b). 



(41) 



Hence, the matching procedure does not lead to C 7 {Mw) and C 7 {Mw) but directly to C| ff (M^) and C^ s {Mw)- 

3. The matching conditions 

The matching for the additional four-fermion operators must only be performed at tree-level. The matching 
conditions are therefore easily derived. An example is given in (pQ). The complete set is 



C Pl (Mw 
Cp (M w 



CfjM- 



w 



CI? (Mw 



Cf 5 (M v 



CfAM w 



CfAM- 



w 



Cf a (M w 
Cf g (M w 



Ct(M w 
Cf 2 (Mw 



1 



1 



R A il2 A jl3 L Lil L Lijo 



\' \'* p e t pe ^ 
& A i22 A j23 l Lil V m m 2 



Gf „ m s m b 2 
T2 KtsKtb ^nU t&n P 



g A il2 A jl3 L LU L Lij m 2 

"1 



1 \/ \/* y.v] pi/ ^ 
g A i22 A j23 L Lii L Lij^T 

1 \/ \/* T-i^f pi/ ^ 



1 



1 



A z31 A j21 L Lil 1 LI, 2 



1 



1 



g-^^^ Lit L£ Jfn 2 



Cf 10 {M w ) = Cf," (M w ) = Cf£(M w ) = 



1 \" \"* p d t pd " 
g A li2^1j3 i fljf 1 Rf-i~^T 



1 \// \//* pdf pd ^ 
- g A H2 A li3 i fljf 1 RZi^T 

dp 



C£(M W ) 
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di 

C< rA M w) = -^2i2^2h T Rje r Rii^r = ^ C< R 3 ( M w) 

Cr 5 (M W ) = Z^12^]l^%^Rli—2- 
° ' n u e 

Cf e {M w ) = -I A ^A^ 3 r^r^-i- = -C%(M W ). (42) 

In Cp^(Mw) we included a term coming from Higgs exchange because it can possibly be large in the case of a large 
tan/3. 

There are some more terms to add to C7, C'$, C7, and Cs, too 



Cy S (Mw) = C7mssm(Mw) 

~ I tf 2k r « Fl ( X <iki>e) 



4m 2 - e 



1 w 



•^fe3i^fe*2i ^R-j£^R£i (43) 



"48 to? fc3i fe2j 



Cg (Mff) — Csmssm(-^w) 
1 1 



4 TO? 



+ ^^T X k3i ) 'k2j r %£ T Re l ( 44 ) 



C7 (Mw) = C7mssm(Mw) 

1 1 



- A ^K a 2^Kl n K nh Y% t Tl l3 [-QuF^x^) + Q e F 2 (x UnSe )] 



e e 



1 1 



+ A ~2~^'ki2^'kj3 T Lz£ r Lij [-QeFl(x ek u ( ) + Q u F 2 (x e 

Up 

4~ m 2 " ^ fc2 ^* k3 ^^ ^ ^ Xdt: P< ^ 



1 '■sd 1 \/ \/* -pd\ r d 

48 m?~ ki2 kj3 LU Uj 



1 1 



^-^Kk2^'jk3^ U Rje^Rli [QdF\{x dk u e ) - Q u F2{xd k u t )] 



2 TO^ 

1 1 1 \ll r df -pd 

2 TO 2 " 2 j3 ^ Mi 



(45) 



(7| ff (Mvy) = C < 8mssm(-^w) 

~ J~T Ka2^%3^in^nbY L ig T e Lej Fi (x u „S e ) 
41 '' t e f 



1 1 



+ 4 ^k*2 ^fe*i3 ^LitFlij ^2 (^efc «<■ ) 
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48 m 2 - fci2 fej ' 3 Li£ L£i 



1 1 



4 m ta 

' _J_\" \"* p^t W 

" 4m 2 ^^fc^ 1 Rjl L Mi 



-2Cp 5 (M, 



H" ) 



(46) 



^. b — ► u|c 



It is convenient to express the branching ratio Br(6 — > S7) through the semi-leptonic decay 6 — > u|c er/ 



T(b -> s-y) 
T(b — > u|c ei/) 



Br(6 -> s 7 ) = v _^- Br exp (b -> u|c eP), (47) 



where we take Br exp (6 — > u|c eP) = 10.5% |54| - This has the advantage that the large bottom mass dependence 
(~ m|) cancels out. In the SM the semi-leptonic decay is mediated by a W^-boson [Fig. || (a)] whereas in the case 
of the MSSM the charged Higgs can take the role of the W. However, the coupling to the leptons is proportional 
to the electron mass and hence it can be safely neglected. Introducing the R-parity breaking terms ( p7| ) -(p9[) offers 
new decay channels depicted in Figs. || (b) and (c) which have to be included in the decay width. Please note the 
following few things: 

• In the MSSM the decay b — > uev is suppressed by the small CKM- matrix element K u b and can therefore be 
neglected. In our case we have to include this decay mode. 

• The absence of lepton generation mixing in the SM forces the anti- neutrino to be v e . This restriction is no 
longer valid in our case. Therefore, we have to sum over the three generations before squaring the amplitude 
because the generation of the neutrino is not detected. 

• Setting the mass of the lepton to zero (which is certainly a valid approximation) the computation in the MSSM 
does not distinguish between the electron and the muon. Here, the two particles involve different couplings. We 
give the results for an outgoing electron. For the muon just change the "1" in the relevant coupling to a "2" . 

The results are then, to leading order and with m 2 /m 2 = 

r(6 ^ s7)= i34 (|C7|2+|(57|2) (48) 

T(b - u\c eu) = { (1 - 8e 2 + 8e 6 - e 8 - 24e 4 log e) x 

[|2,4 + C 2 | 2 + |5 2 | 2 ] +|Bi| 2 + |C 1 | 2 }, (49) 



where 



A = 2V2G F K 23 (50) 

o _ \ " ^»jl^mn3 r e pef K r— 1 2 (5Pl 

i=l 
3 w \/» 

r* _ A i3k A lmn - r d pdf /rr,\ 

i—1 dp 

e = — . (53) 
mb 
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IV. RESULTS 



The formulas of the previous section are far too complicated to be treated "by hand" but it is no problem to feed 
them to a computer. The 7-matrix is independent of the parameters of supersymmetry. With the help of Mathematica 
p5| it is possible to diagonalize it and find the influence of the QCD effects. 



A. General results 



Because of the enlarged operator basis the expression for C7 changes. In general, the solution of the RGE for the 
Wilson-coefficients is given by 



a s {M w ) 
a s {p) 



TP 
Hi 



where V diagonalizes 7 T 



Id = V'^V. 



Po = 23/3 is the one- loop beta- function and 7d is the vector containing the eigenvalues of 7. In our case 

-a, i-ii) -16 -16 -16 -16 



(-16 
-16 
4 



-16 
4 



-16 

4 



(54) 

(55) 
(56) 



4 28 28 

3 3 



-VI 

3 



f 2.233 



2.233 6.266 6.266 -13.791 -13.791 -6.486 -6.486). 

To have an idea which coefficients are relevant we perform a numerical ana-lysis with a s (Mz) — 0.121 and /i 
m& = 4.2 GeV. The coefficients C^ s (m,i)) and C^ s (m b ) are then 



Cf(m b ) = -0.351C^(Mw) - 0.198C%{M W ) - 0.198(7^ (M w ) 
-0.l7&Cf B (M w ) + 0M5Cf(M w ) + 0.093Cf(M w ) 



Cf(m b ) = OMOCf^Mw) + 0M0Cf 2 (M w ) - 0.198Cf,f (M w ) 
-0A98CfjM w ) - 0.178C C P S JM W ) + 0.381C^(M W ) 
+0.381C^(Mtv) - 0.213C^(Afw) + 0M5Cf(M w ) 

+0.093Cf(M w ), (57) 

It is clear that the four-fermion operators including a left-handed s-quark contribute to CV/g whereas the ones with 
a right-handed s-quark contribute to C7 /§ . The numbers multiplying the different Wilson coefficients are all of the 
same size, hence there is a priory no term which can be neglected. 



B. Specific results for R- parity violation 

It is obvious that the Wilson coefficients depend in a very complicated way on the parameters of our supersymmetric 
model. Changes of A, A', or A" not only affect the result in a direct way but also in an indirect fashion through an 
altered mass spectrum and different mixing matrices. Therefore it is very hard to make general statements on the 
behaviour of the branching ratio. 

As mentioned, it is not our aim to perform a high precision analysis of the parameter space [pl^-jsl)) but to explore the 
influence of the R-parity violating couplings on b — > sj with a reasonable accuracy. We know that our results for the 
branching ratio are only valid at the 25% level because we do a leading-log approximation with large scale uncertainty 
p9[ . However, we expect the impact of the R-parity breaking terms not to change much when including the next- 
to-lcading corrections. This means that the shape of the curves remain more or less the same whereas the offset 
where our curves start (i.e., no R-parity breaking) may change significantly when calculating the next-to- leading-log 
approximation. Solving the RGE's was performed in the following way: 
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• Solve the equations for the gauge couplings. The boundary conditions are the physical values at Mz- The gauge 
couplings will meet at Mx ~ 2 • 10 16 GeV. Choose a common value M(Mx) for the gaugino masses at that high 
scale and solve the RGE's for M1-M3. 

• Set tan/3 and X(Mz), X'(Mz) and \"{Mz) to the desired value and use them together with the quark and 
lepton masses as inputs at Mz for the Yukawa couplings. Let these couplings run to Mx- 

• Choose A, fj,2i(Mx) and trial values for m and fi(Mx) to complete the boundary conditions at Mx- The whole 
set of RGE's is then run down to Mz- We choose /J.2i(Mx) = 0. This is an approximation because /i 2 i and fi2i 
will not vanish at Mz- However, in our examples their values are so small compared to /1 and jl that we can 
neglect them avoiding a mixing between H\ and Li. 

• The minimum of the one-loop effective Higgs potential Vq + AV will in general not be at (wi,«2) which makes 
it necessary to adjust /i(Mx) and m in a clever way. We chose the Newton method to converge to the desired 
position of the minimum. Following [ p"8| we evaluate the minimum of the potential at some average mass scale 
to avoid large logarithms and therefore get a more reliable result. 

The next step is the numerical diagonalization of the mass matrices to find the masses of the physical particles and 
the relevant mixing matrices. 

Changing the value of the R-parity violating couplings makes it necessary to continuously adjust the values of uiq 
and fi(Mx) which alters the mass spectrum of the particles at Mz- All the following examples correspond to mass 
spectra within the current bounds for the masses of the supersymmetric particles. We encountered two critical 
situations, namely too small masses for the lightest selectron and/or the lightest Higgs boson. An idea would then 
be to constrain the bounds on the R-parity violating couplings through the requirement of a phenomenologically 
realistic mass spectrum independent of the value for the branching ratio Br(6 — + S7). However, in general, a realistic 
mass spectrum does not restrict the R-parity breaking parameters substantially. Moreover, the mass spectrum highly 
depends on the values of tan/3, A and M(Mx)- Hence, such bounds would be strongly model dependent. 
As reference model we chose tan/3 = 5, A = and M(Mx) = 300 GeV. With vanishing R-parity violating couplings 
this leads to squark masses of 600 — 800 GeV, sleptons of 260 — 330 GeV, a lightest neutralino of 120 GeV and a 
lightest Higgs of about 100 GeV. Fig. || shows the behaviour of Br(6 — > sj) in the neighbourhood of our reference 
model with, in addition, A' 132 = A' 122 =0.1. Interestingly, the value for the branching ratio is rather stable under a 
change of tan (3 and A. 

As mentioned before, it is very difficult to isolate generic features of the different models. Let us make some comments: 

• At least two of the A's must be non-zero to have an influence on the result. There are two exceptions: A" 2 3 and 
A 223 alone will give a contribution due to the anti-symmetry of A". However, their impact on the branching 
ratio is so small that no reasonable bounds can be found. 

• As all the effective couplings depend on the inverse mass-squared of a heavy particle, it is clear that the influence 
of the new physics is bigger in models with a lower mass spectrum. To see the effect of smaller masses compare 
Figs, [n^ (a) and (b) where for the latter M(Mx) = 100 GeV is taken which results in, for instance, squark 
masses of about 250 GeV. Note, that our reference model leads to quite high masses. 

• The squarks are always heavier than the sleptons. As a consequence, the influence of non-vanishing A" on 
Br(£> — > sj) is much smaller than for non-vanishing A'. 

• The simplest non trivial situation consists of a single non- vanishing pair of R-parity violating couplings. We 
encounter the following scenarios: 

— The largest effect results from the pair (A' 122 , A' 132 ) (see Fig. [ll]). It is not reasonable to extract a bound 
for the product of these parameters because (through the RGE's) the dependence on these couplings is 
much more complex. However, we draw the conclusion that significant effects result if both couplings 



exceed 0.1-0.2. Fig. 11 also shows that the branching ratio may strongly depend on the relative sign of the 
couplings. This happens if the new contributions are able to diminish the value of C7. 

A general feature of our model is that the influence of a pair of non- vanishing R-parity violating couplings 
on the branching ratio starts being significant if (A', A') ~ 0.2 — 0.4. For (A", A"), in most cases, the 
requirement of non-diverging RGE's for these couplings gives more stringent bounds. There is one more 
aspect: Some of the contributing pairs have more than one distinct index. They appear only because of 
the mass differences of the different squark or slepton generations. This can be seen clearly in formulas 
!2h: If the masses were independent of their index the T's would combine to an identity matrix leaving 
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only those pairs with an identical sfermion index. The bounds on these pairs are much less stringent. This 
is exactly what happens in the case of X'{ 2 3 and X'223 when being the only non-vanishing coupling. This is 
another reason why the results for this situation are not stringent. 

— A special situation is depicted in Fig. The branching ratio seems to "explode" at A^ 12 = A^ « 
0.3. This is due to the mentioned decreasing mass squared of the lightest selectron, which appears in 
some denominators in the matching conditions. Hence, the specific position of the peak is highly model 
dependent. 

• If several pairs of couplings are non- vanishing the picture can get more complex, as Fig. [[3] shows. There, we 
combined the decreasing effect with the peak due to the small selectron mass. 

• It is difficult to compare our results with existing constraints for the R-parity breaking couplings because the 
results are extremely model dependent. However, in the case of A' our bounds are highly competitive. For 
comparison we put the current bounds pl-pll in the captions of the respective figures. 



V. CONCLUSIONS 



To examine the influence of R-parity breaking on b — > S7 one has to enlarge the operator basis substantially. At 
the leading-log level it consists of 28 operators, neglecting Higgs-lepton mixing avoiding this way the introduction of 
scalar operators. The corresponding 7-matrix can be found with the help of previously known results and diagonalized 
numerically. The matching conditions of the magnetic penguins Oj and 0$ get new contributions. Their counterparts 
of opposite chirality, O7 and Os, also have to be considered. If one uses the semi-leptonic decay b —> u\c ev to cancel 
the large bottom mass dependence new contributions to this decay must be included. 

R-parity breaking definitely has influence on the branching ratio of 6 — > 57. However, its impact is highly model 
dependent because the (unknown) supersymmetric masses are mostly responsible for the size of the new contributions. 
In a cautious model the new couplings are able to change significantly the result if they are of order lO^ 1 . Moreover, 
45 new (complex) Yukawa couplings offer an infinite number of possible scenarios. The simplest cases involve only 
one or two couplings present but also in these situations completely different evolutions of the branching ratio are 
possible. To give more concrete results we definitely need more informations on the parameters of our supersymmetric 
model. 
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APPENDIX A: MIXING MATRICES AND INTERACTION LAGRANGIAN 



1. mixing matrices 




In this first appendix we present the mass mixing matrices for the relevant particles. They are needed for two 
reason: First, their eigenvalues correspond to the physical masses of the particles and second, the unitary matrices 
that diagonalize the mass matrices rotate the fields to their (physical) mass eigenstates. 



Charginos 



The charginos X1/2 are a mixture of charged gauginos X^ and Higgsinos h 1 and h\. Defining 

* + =m ^-(T) < ai) 
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the mass terms are then 



Ct = --{i> +T X T ^~ + ^- T X^+) + h.c. 



where 



X 



-M 2 92V2 



92V1 -H 

The two-component charginos \f (i — 1, 2) and the four-component charginos \ f/ 2 are then defined as 



xt = I — I xt 



xi 

Xi J V X2 

where the unitary matrices U and V diagonalize X: 

m£ = u*xv- Y = vx^u T 

Ct then becomes 

r ch — — M ch v cll v ch A/f ch -i/Ch ch 
J -m — N1 D\\X\ Xl ~ 1V1 D22X2 X2 ■ 



(A2) 



(A3) 



(A4) 



(A5) 
(A6) 



U and V can be found by observing that 



M£ h2 = vi T iy _1 = u*xx T u*-\ 



(A7) 



They are not fixed completely by these conditions. The freedom can be used to arrange the elements of M|, h to be 
positive: If the i th eigenvalue of M|, h is negative simply multiply the i th row of V with —1. 



Neutralinos 

The neutralinos are linear combinations of the gauginos A' and A 3 and the neutral Higgsinos h\ and h®- If we define 




the neutralino mass term reads 



where 



V> 0T 1>° + h.c, 



Y = 



Two- and four-component neutralinos must be defined as 





-Mi 





givi 


3l«2 
72 


\ 







-M 2 


92V1 

V2 


92 V2 

V2 






9l fl 


S2«l 











72 






V 




92f 2 
72 


M 





/ 



2=1, 



(A8) 



(A9) 



(A10) 



(All) 
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To diagonalize the mass matrix N must obey 

N D = N*YN~ 1 , (A12) 

where Nd is a diagonal matrix. N can be found using the property 

N% = NY^YN' 1 . (A13) 

The eigenvalues and eigenvectors are found numerically. Possible negative entries in Njj are turned positive by 
multiplying the corresponding row of TV by a factor of i. 



Quarks and Leptons 

The situation in the quark and lepton sector is in almost complete analogy to the standard model. The quarks and 
leptons get their masses from the Yukawa potential when the Higgs bosons acquire a vacuum expectation value. We 
define the mass eigenstates by 

•J? = ut 3 u Lj «.&> = un URj 

d£?> = D^d Li dft> = D%d Rj (A14) 
(™) _ F L e (m) -E R ev 

The mixing matrices must satisfy 

D R X dT D L^ = \d D = diag(^f) i = l,...,3 

u R X uTtjL] = A u = diag ^I^^ ( A15 ) 
E R\eT E L] = ^ = diag 



v\ = V 2 cos p 

92 

v 2 = V2— sin/3. (A16) 

92 

As one can see, the eigenvalues of X u and X d are fixed by the quark masses and the minimum of the Higgs potential. 
In the SM the only effect of the mixing which can be seen is the CKM-matrix K — U L D L ^ appearing in the flavour 
changing charged currents. Therefore it is possible and convenient to set 

D L = D R = U R = 1 {=> U L = K) (A17) 

(To be more precise, one chooses X d and A e to be diagonal and A" = K T di&g (m U i / v 2 ) ■) Although in our theory the 
mixing matrices appear in all kinds of combinations we adopt this convention here, emphasising that it is a choice 
made just for convenience. It is possible that one day an underlying theory fixes the values of A" and X d at some 
(high) scale. 

Please note that in the text we neglect the superscript m for the mass eigenstates. 



Squarks and Sleptons 

If supersymmetry were not broken squarks and sleptons would be rotated to their mass basis with the help of the 
same matrices as their fermionic partners. But since this situation is not realistic we need to introduce a further 
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set of unitary rotation matrices. The notation must be set up carefully because the mass eigenstates of squarks and 
sleptons are linear combinations of the partners of left- and right-handed partners of the corresponding fcrmions. We 
define the 6 x 3-matrices T in the following way: 



(A18) 



S(m) = ( V e Ipe \ { e L 



To diagonalize the mass terms the mixing matrices have to satisfy 

r ^ u r«t = m 2 uD 

rV d r dt = m s 2 dD 
r e m s 2 c r £ t = m 2 cD 

T»mlr^=ml u , (A19) 
where the matrices on the RHS are diagonal containing the masses squared of the phy-sical particles. 
The mass matrices are (with the exception of the sneutrino) of the form ^ _gt (7 ) ' w ^ ere ^' ^ an< ^ ^ are 3x3- 
matrices. For the different fields they are (choosing fi real) 
• up-squarks: 

A = m 2 + w 2 A"A Mt + Mf cos 2/3 ( 1 - § sin 2 W ] 1 3 



2 3 



B =(h u + ficot(3X u )v 2 

C = m 2 u + v%\ u ^\ u + \m% cos 2/3 sin 2 6 W 1 3 (A20) 



• down-squarks: 



A = m 2 Q + v 2 1 X d X d ^ + Mf cos 2/3 ^-1 + 1 sin 2 6>v^ la 



B = (h d + f i tan f3X d ) i>i 



sclcctrons 



The sncutrinos have 



C = m 2 D + vjX d ^X d - \m 2 cos 2/3 sin 2 %1 3 (A21) 



A = m\ + w 2 A e A et + M| cos 2/3 + sin 2 6 W J 1 3 
B = (/i e + ^itan/3A e )wi 

C = m| + w 2 A et A e - M| cos 2/3 sin 2 (9 W 1 3 (A22) 



m\ v =mj + ^- cos 2/31 3 (A23) 
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Higgses 



The Higgs sector consists of two SU(2) doublets hi and h 2 - The real and the imaginary part of the neutral 
components mix via the matrix 



[i 2 + mjj + Si^2.(3v( - vi) 

1 2 i 2 



2, 2 

~ _ 9±_92_ ViV2 



M 2 + m 2 H 2 



(A24) 



and 



M 2 + m 2 H2 



(A25) 



respectively. The charged components give rise to a mass matrix of the form 



^ +m 2 H +9 1 +9 3 . {v 2_ v 2 )+ 9 3 v a . 

A, 



-fl + -fviv 2 

~2 1 2 



h -i- « in<it„ // 2 -I-tt7 2 9l+32 f,,2 2\ 1 g 2 ^i 

-H+—V1V2 fl +m, H2 j — \ v i — v 2)~' T~ 



(A26) 



All the above expressions are valid at tree-level. The vacuum expectation values v± and u 2 are found by minimizing 
the e ffecti ve Higgs potential. If one inserts the gained formulas into (A24) - ( A26| ) one of the eigenvalues of ( A25| ) 
and ( A26 ) vanishes, indicating the eaten fields of the Higgs mechanism that takes place. 



2. Interaction Lagrangian 

For the evaluation of the matching conditions we need certain parts of the interaction Lagrangian. In addition to 
equations (p?]) - (p9|) these are 



S 'quark- Quark- Chargino 



C- qq ^ = ufdt [Af je P L + Bf je P R ] xf c + u\xf c A%P R + B%P L 



+d j u i [A u ijt P L + B? jt P R ] xt + d\xf 



A?],P R + B?},Pr 



J ■■ 



(A27) 



where 



A" 



d r u U..r[* 



{\ d D T 



Bt 3l = {K^\lT R %V l2 ~g2T u L \ 3 Vn 



A? jt = (\ u D KT d L %V e * 2 



m t = (K\ d D r R % j u e2 - 92 (KT a I ;) ij u a 



Pl/r = 

and xf 1 ° denotes the charge conjugated field. 



(A28) 



Cqq X 



S quark- Quark- Neutralino 
-djU [C d e P L + D d l3l P R ] X °i ~ d{xj 
-tijWt [C% t P L + D% t P R ] X ° e 



C%Pn + DfJ e P L 



C^Pr 



D$.Pt 



Uj, 



(A29) 
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where 



Ct Jt = (^ d L %N £ * 3 - V2giQ d T%N a 



D'i 



U = frM^Nn + ±=T%{(2Q d + l) gi N n - g*N a ) 



D%t = (^ U J)ijNu + -j={KT^M{2Q d + l)gxNa+g2N l2 ). 



(A30) 



S quark- Quark- Glumo 
C iqg = -V2g s T A ^ [u\ a gA [(TIK%P L - T RlJ P R ] u jP 
+ dig! [T d Uj P L - r d mj P R ] djp 



(KT^) tJ P R -T R \ 3 P L 



fjA 



1 Lij^R - 1 m-i^L 



9 A 



(A31) 



Gluino- Gluino-Gluon 
■ \9sJ ABC 1Ta1,9bA» c 



(A32) 



Note: There is a symmetry factor of two in the Feynman rule for this vertex. 



APPENDIX B: RENORMALIZATION GROUP EQUATIONS 

We present the full set of RGE's for all the parameters of the MSSM including R-parity breaking terms. Our results 
are in complete agreement with |6C]]26|] , although we don't restrict ourselves to couplings of the third gen eration. All 
the formulas can be derived from the expressions for the most general form of a softly broken SUSY |6lj |63| . Let us 
begin with the parameters of the superpotential W [t = In /z) . 



16^ 2 - M = Tr (3A d A dt + 3A u A ut + A e A et ) p - [gj + 3 5 |] n 
+ [3A^ m A^ m + A| m A*,g m ] fi2i 



(Bl) 



^k 2 —^ - [(A e A et ),y + A« m A^ OT + 3A^ m A^ m ] fj, 2 j 



+3Tr(A"A Mt )M 2l - [g{ + igl] M2i 



(B2) 



dt 13 



[(X u X u ' r )ik + (A d A dt ) jfe + X' e * km X' Um ] \%j 
+ [2(X U X v )kj + X'li m X" em ] X? k 
+3Tr(A"A ut )A^ 



13 , „ n 16 , 

g-ffi + 3g 2 2 + y g\ 



1671-2 ^t^' - [(A"A" 1 ') i fc + (X d X df ) lk + X'[ km X' Um ] Xi 3 



(B3) 
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+ [2(A dt X d ) kj + 2\' e * mk \' emj + 2\" km \" jm ] \k 



[3Tr(A d A dt )+Tr(A e A et )] A," 
[3A^ m A' 



/,-/ i + A £m A ktm 

16 



a; 



kij 



\a\ + 3s! + y al 



[(A e A et ) jfc + \l em \um + 3Aj £ * £m A- £m ] A^ 

+ [ 2 (^ et A e )fcj + Kmk^lrnj] Kk 

+ [3Tr(A d A dt )+Tr(A e A et )] A?. 

+ 2 [3^m^fc£m + ^Im^Wm] 



167T -A iifc 



[(A L A et )i„ + A* £m A^ m + 3A^* £m A- £m ] A„jfc 

+ [(A e A e t)j„ + A*£ TO Aj^ m + 3A^ TO A^ m ] A iri fe 
+ [2(A L ^A e )„fc + A| mn A£ m fe] Ay n 
+2 [3A^ TO A i£m 4- A £m Ai£ m ] A^ fe 
— 2 [3A^A^- £m + Af^A^m] A? fe 

- [3gl + 3#f] Ayfe 



16^ 2 4-A' 



dt 



^ijk 



[(A e A L ^) in + A* £m Ai£ m + 3A^ £m A- £m ] X' n j k 
+ [(A u A ut )., n + (A d A dt )j- n + X' e * nm X' ejm \ \' ink 
+ [2(A A + 2A^ n A £mfc + 2A"* m A" fcm ] A- Jn 



[3AfcjA^ m 4 

+ 3 5 2 



16 



"33 



'Hjk 



df yfc 



[2(A"^A") n , + A" £m A" £m ] A" j7c 
+ [2(A dt A d ), 



^ A lnm A lmj ' ^ A lnm A tjm\ A ink 



+ [2(A dt A d )„ fe + 2Xf nm 

A imk ^ A lnm A lkm\ A iji 

"4 



A" 

/s ijk 



The parameters of V^ ft obey the following RGE's: 

d 



dr 



Tr (3A d A dt + 3A u A ut + A e A et ) fl 

+ \^ A im A 'ilm + Km A iem] M2i 

+2Tr (3/i d A dt + 3/i"A ut + h e A et ) p 

+2 [3^ m A- £m + frf m A* £m ] [i2% 

- [al + 3g$\ M + [2M ig l + QM 2 gl] 



167T -fo 



[(A e A e % + X*e m Xu m + 3A^ ro A^ ro ] £ 2j + 3Tr(A"A Mt )M2, 
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+2 [(/i e A et )y + C'u m X* em + 3C' iem \'* gm ] /j,2j 
+6Tv(h u X u ^ 2i + 2 [3\tfiu m + K* m C Um ] fi 
- [gl + 3.9a] fci + [2M l9 l + QM 2 gl] fi 2l 

16tt 2 ^^ - 5(A"A Mt /i u ) y + (X d X d ^h u ) ij + 4(h u X^X u ) l3 

+2(h d X^X u ) t3 + 3Tr(A u A ut )/i% + 6Tr(/i"A ut )A£ 

\' h u -I- \"* \" h u -l-9\'* C 1 \ u 
+ 2 ^emn C jmnKe + "g^l (2-^lA^ ~ Kj) 

+3g 2 (2M 2 A« - h%) + y.g 3 2 (2M 3 A^ - h%) 

^jh% = 5(A d A dt /i% + (X u X u ^h d )ij + 4(h d X d ^ X d )ij 

+2(^A ut A% + Tr(3A d A dt + A e A et )/i^- 
+Tr(6/i d A dt + 2h e X^)X d tJ 
+ [^em^'kij + 2/i^Aj^] X* Mm 
+3 [^/m^ij + 2h d m X' ki j] X' k \ m 
+ [^fejAft m + 2\ k jC' eim \ X'l km 
+2 [^ife^mj + 2A ife C^ TOJ -] A 
+2 [^ifcA^m + 2A ifc C^ m ] A 
7 

+ 9 

+ -gl{2M z X d 3 -h d 3 ) 

16tt 2 ^^. = 5(A e A et /i%' + 4(/i e A et A%- 

+Tr(3A d A d t + A e A e t)^. + Tr(6/i d A d t + 2h e X e ^)X e l0 

+ [A? m Cfeij + 2/l| TO Afejj + h%jXu m + 2X% j Cum] Kim 

+3 [A^ m Cfcy + 2h d m Xkij + h%jX' i(m + 2X e k jC' ilm \ X klm 
+3. 9 2 (2M 1 A^- - ft? .) + 3. 9 2 (2M 2 A^ - ft?.) 

167T 2 — Cjjfe = [(A e A e ^) in + A* £m Ai£ m + 3A^ TO A-£ m ] C„jfc 



1 //* 

Vfcm 

5 2 (2MxA^ - + 3. 9 2 (2M 2 Af :) . - 



+ [(A L A et )j„ + A* £jn Aj^ m + 3A^ m A^ m ] C in k 

+ [3A^A^ m + Af^A^ TO ] hj k - [3A^*„A^ m + Af^A^ m ] h 

+ [2(A e ^A e )„fe + X* tmn Xi m k\ Cij n 

+ \2{h e X e ^)i n + 2X* nim Cu m + 6X'* lm C' iem ] X n jk 

+ [2(h c X eji )j n + 2X* ntm Cji m + GX'* em C' jem \ Xi n k 

+ [^^tm^Um + 2A|^C^ TO ] X C j k 

~~ [S^emCjem + 2A^Cj£ m ] X e ik 

+ [4(X C ' f h C ) nk + 2X* lmn Cl m k\ Xij n 

+3gl(2MiXij k — Cijk) + 3<7 2 (2M 2 Ajjfc — C^k) 
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16n 2 — C' ljk — [(A e A et )i„ + X* nlm Xn m + 3A^* £m A^ m ] C' njk 



- [(X u X u ^)j n + (A d A dt )j„ + A^ m A^ m ] C,' nfe 

( ink 

' [3A^A^ m + A|^A^ TO ] hj k 

- \2{h e X e ^)i TL + 2X* nlm Cn m + 6X'* em C' iem ] X' n j k 

- [2(h u X^) jn + 2{h d X d ^) jn + 2X'l nm C' ejm ] X' mk 

- \A(X d ^h d ) nk + 4:X' e * mn C' emk + 4A"* m C" fcm ] X' ijn 
' [6A^C-£ m + 2X^Ca m ] X d k 



-g\{2M x X\ ]k - C[ ]k ) + 3g 2 2 (2M 2 X' ijk - C' ijk ) 



+ l ±g 2 (2M 3 X' ijk -C' ijk ) (B15) 
167T 2 — Cy- fe = [2(A ut A")„, + X"* hn X" em ] C'n jk 



+ [2(X d ^ X d ) n j + 2X'l mn X' lm j + 2X"* m X"j m ] C" nk 
+ [2(A dt A d ) nfc + 2X' e * mn X' emk + 2A"* m A" fcm ] C" jn 
+ [4(A ut /i u ) ni + 2X"*e m C" em ] X'n jk 
+ [4(A d ^/i d )„j + ^X'l mn C'i m j + 4A"* m C^ m ] X" nk 
+ [4(A dt /i d )„ fe + 4:X' e * mn C' emk + 4A"* m C" fcm ] X'- jn 

+ \g 2 {2M 1 X'( jk - C'r k ) + 8g 2 (2M 3 X? jk - C'l jk ) (B16) 



l^jml^ = [2(A»tA") w + A- m A^ m ] m 2 UkUj 



+ [2(X u 'X u )jk + X"l m X' klm \ m 2 JtUk 

+ 4A ^ A mi m Q t Q m + q:A ji A U m H 2 H 2 + 4A j£m A mm m D,fl„ 

+4(/ l "V) + 2C$? ro C£ m - -^-g 2 M 2 5ij - -glMl5 t3 



167T 2 — m|, i£ ) 3 — [2(A dt A d )jfe + 2A^ J -A £mfe + 2A"* m A" fcm ] m 2 D . Dk 



+ [2(X dJf X d ) ki + 2X'l mk X' tmi + 2X" km X" im ] m 2 DkD . 
+<±Xp d u m 2 HiHi + 4Xf e X d mi m 2 QeQm + 4X'; nj X' mnl m 2 LeLm 
+4X'* e jX' nmi mQ e Q m + 4:X^* ej X" mi m 2 DeDm + 4A"* J A^ ni m 2 w 
+ AX d n *X' kni ml k 

11 I, 

32 
~3 



+ [(A"A»t) ifc + (A d A d t) tfe + X>; km X> m ] m 2 QkQj 
+2A"j X™ e m 2 H2H2 + 2X^ 1 \ j Xy t m 2 , mUe + 2Xfj X d l m 2 HlHl 
+2Xj tj X d e m 2 DmDe + 2X'* ljn X' ein m 2 LmLe + 2X'* jm X' nil m 2 DmDe 



(B17) 



g 2 M 2 5 l3 - -g 2 3 M 2 S t3 (B18) 
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+2(/ l »ft»t ) . . + 2 (h d h% j + 2C'; jm C' Um 



2 S2 
- -ffiMl&n - GglM^ - — g 2 M 2 5 l3 (B19) 



16^^!^. = [2(A fi tA e ), fe + \* emj Xemk] m% iEk 
+ [2(A et A e ) fe4 + X} mk X tm i\ m EkEj 
+4Xf k X e ki m 2 HlHl + 4\f m X e H m 2 LmLe + 4X* nmj X nU m 2 LmLe 



+^KkjKi m L k H 1 +4A, e *A nfcl m 2 



H\L k 



+A{h^h e ) ij + 2Q mj C emi - 8g 2 M 2 S tJ (B20) 



2 



16^ 2 ^m| iLj = [(A e A e t) 4fc + A^ m A« m + 3A- m A^ m ] m LkL . 

+ [(A e A et ) fe:) - + X* j£m Xke m + 3X'* em X' kem \ m 2 L . Lk 
+ \3XimKim + A|mA^ TO ] m 2 HlL . 
+ [3A^ m A^ m + A^ m A| TO ] rn 2 L . Hi 

+2X c l } j X l i k m 2 HiHi + 2X e ^X e it m 2 EmEe + 2X* mjn Xu n rn 2 LmLe 
+2X^ : j m X n um 2 EmEt + 6X'* mn X' Un mQ m Q e + 6X'* nm X' ine m, 2 DmD( 

+2^kjnKn m L k H 1 + 2A^* Afcin^ 2 ^^ 

+2(h e h eji )i j + 2C* em Cu m + §C'* lm C' iS , m 

-2g\Ml^ - 6glM 2 S tJ (B21) 

l^ 2 jm 2 HiHi = Tr(6A d A d t + 2X e X^)m 2 HlHl 

+ [3A-£ m A^ m + X* Um X e em \ m 2 L . Hl 

+ [SXgmKem + A|mA«m] m 'H 1 L i 
+6^im^en m Q m Qi + ^mn\l m 'b m D t 
+ ^^nmKn m L m L t + 2A^„A^m|; m£{ 

+Tr(6/i'V t + 2/i e /i et ) 

-2 3 2 M 2 - 6. 9 2 M 2 (B22) 



+Tr(3A d A dt + A e A et )m| iHi 
+ [3A^A^ m + A^Ai£ m ] m 2 



+ [3A^A' fc£m + A^A Mm ] m| iLjc 
+6A^A^m|, mr , < + 6X^ n X' an mQ m Q t 
+ 2 K*m X nUm 2 EmE£ + 2X e * m X Un m 2 LmLil 

+6hf^ n C' iem + 2h e ^ n Cu m (B23) 



lQK 2 jm 2 H2H2 =6Tr(X u X^)m 2 H2H2 



+®Kirn X en m Q m Ql + ^KnnKit. m U m U t 

+6Tr(/i"/i ut ) - 2 3 2 Af 2 - &g 2 2 Ml (B24) 



2G 



Gauge couplings and gaugino masses run in the following way 

16^ 2 ^32 = 9\ 167r a ^M a = 2g\M 2 (B25) 



APPENDIX C: THE 7-MATRIX 



The Q CD-mix ing of our new operator basis leads to a 28 x 28-matrix. How its elements are deduced is explained 
in section pll B 2 . Fortunately, many entries vanish giving us a chance to der ive the eigenvalues and -vectors with the 
help of Mathematica. We split j cS in the three block mentioned in section IIIB2. The four-fermion operators give 
the block (we have included the ordering of the operators in the first row) 

/O1O2 O3 4 5 Oe, Pi P 2 P3 Pi P5 Pe P7 Ps P9 P10 P11 P12 Ri R2 R 3 Ri R5 Re\ 
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The 4 x 4-block of the magnetic penguins looks like 

/ o 7 o s o 7 



32 
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_32 
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28 o 
3 u 
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_ 32 
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(CI) 



(C2) 
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The four last columns that mix four-fermion operators with magnetic penguins are given as rows. They are 
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(C3) 



We emphasize that the matrix depicted here is 7 . In the HV-scheme it should coincide with the uncorrected 7. We 
have checked this explicitly for all the entries. 

APPENDIX D: DEFINITION OF THE FUNCTIONS Fi - F 4 

These functions appear in all calculations of diagrams like those of Fig. |6[ 

= -10/1 ~ ^4 [2 + 3x - 6x 2 +x 3 + 6x\ogx] 
12(1 — xp 

F ^ x ) = TTTTT- — vT [1 - 6ar + 3a; 2 + 2x 3 - 6x 2 log x] 
12(1 — xp 

Fz{x) = 2(1 - xf + 4X _ X " - 2l °Z x 1 

F4x) = —l—-[l-x 2 + 2x\ogx] (Dl) 
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FIG. 1. From the full to the effective theory 




FIG. 2. Diagrams contributing to the 1-loop mixing of the four-fermion operators. 




b s 



FIG. 3. Diagram that gives finite contributions. 




FIG. 4. Two-loop diagrams needed for calculating the 7-matrix. The wavy line can be a photon or a gluon. These diagrams 
have no closed fermion loop. 
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FIG. 5. Two-loop diagrams needed for calculating the 7-matrix. The wavy line can be a photon or a gluon. These diagrams 
have a closed fermion loop. 




b W s b H s b u s b d s b d s 



(a) (b) (c) (d) (e) 

FIG. 6. Diagrams that contribute to the matching of O7/8 and 0?/g. The outgoing photon/gluon is attached at every 
possible position. 




(a) (c) (b) 

FIG. 8. Contributions to b — * u\c ev. The arrows indicate the particle flow. 
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FIG. 9. Behaviour of Br(& — > s'y) in the neighbourhood of our reference model with tan/3 = 5, A = 0, p,2i = 0, M = 300 
GeV and, in addition, X' 132 = A'122 = 0.1. The horizontal lines show the current experimental bounds. 




(a) (b) 



FIG. 10. Values for Br(6 — > S7) as a function of A212 = ^213- The other parameters are: tan/3 = 5, A = 0, p,2i = 0, and (a): 
M = 300 GeV, (b): M = 100 GeV. The current bounds are: A212 < 1-23, A 2 ' 13 < 1.23 Q, IA212 • A 2 ' 13 | < 0.006 10. 




The current bounds are: 
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0.1 0.2 

y — v 

A 312 — A 313 

FIG. 12. Values for Br(b — > sj) as a function of A 312 = A 313 within our reference model. The current bounds are: 
A312 < 0.11 10 ^q cV , A 313 < 0.11 10 ^q cV JlH , |A 312 ■ A313I < 0.01 for squarks of 100 GeV [Q. The peak results from a small 
lightest selectron mass. 
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A 122 — A 132 — A 321 — A 331 — ZA 312 _ A 313 

FIG. 13. Values for Br(6 — > S7) as a function of — A' 132 = A122 = — A 331 = A 32 i = 2A 312 = A 313 within our reference model. 
Additional current bounds are: A331 < 0.45 10 o^ cV 1 A^ 2 i < 0.52 10 ^ eV [ pr[ . 
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